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ABSTRACT 

o" 

■ We consider hot, two-temperature, viscous accretion onto a rotating, unmagnetized 

CN . neutron star. We assume Coulomb coupling between the protons and electrons, and 

\ free-free cooling from the electrons. We show that the accretion flow has an extended 

^ ' settling region which can be described by means of two analytical self-similar solutions: 



a two-temperature solution which is valid in an inner zone, r < 10^'^, where r is the 
radius in Schwarzchild units; and a one-temperature solution which is valid in an outer 
CN ■ zone, r > 10^'^. In both zones the density varies as /) oc and the angular velocity 

as VL on r~^/^. We solve the flow equations numerically and confirm that the analytical 



> 



^ . solutions are accurate 

o 



Except for the radial velocity, all gas properties in the self-similar settling zone, 
^) \ such as density, angular velocity, temperature, luminosity, angular momentum flux, 

' are independent of the mass accretion rate; these quantities do depend sensitively on 



o 
o 
o 



the spin of the neutron star. The angular momentum flux is outward under most 
PLh' conditions; therefore, the central star is nearly always spun-down. The luminosity 

Q ■ of the settling zone arises from the rotational energy that is released as the star is 

braked by viscosity, and the contribution from gravity is small; hence the radiative 
. efficiency, rj = Laccl^(? ■, is arbitrarily large at low M. For reasonable values of the 

gas adiabatic index 7, the Bernoulli parameter is negative; therefore, in the absence of 
dynamically important magnetic fields, a strong outflow or wind is not expected. The 
' flow is convectively stable, but may be thermally unstable. The described solution is not 

advection-dominated; however, when the spin of the star is small enough, it transforms 
smoothly to an advection-dominated branch of solution. 

Subject headings: accretion, accretion disks — stars: neutron 



X 



1. Introduction 

At mass accretion rates less than a few per cent of the Eddington rate, black holes (BHs) are 
believed to accrete via a hot, two-temperature, radiatively inefficient, quasi-spherical, advection- 
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dominated accretion flow, or ADAF (Ichimaru 1977; Narayan & Yi 1994, 1995a,b; Abramowicz et 

al. 1995; Chen et al. 1995). The physical properties of ADAFs around BHs have been investigated 
by a number of authors, and detailed spectral models have been applied to observations of BH 
candidates (see Narayan, Mahadevan & Quataert 1998 and Kato, Fukue & Mineshige 1998 for 
reviews). 

At low mass accretion rates, accretion onto a neutron star (NS) is also expected to occur via 
a hot, two-temperature flow (Narayan k. Yi 1995b; Yi et al. 1996), but the properties of such flows 
have not been investigated. NS flows are expected to differ from BH ADAFs in several respects, 
(i) Whereas in the case of a BH the accreting material flows freely and supersonically through the 
absorbing boundary at the event horizon, in the case of a NS the radial velocity of the material 
must decelerate to zero, (ii) The accreting material is expected to apply a spin-up or spin-down 
torque on the NS. Popham k. Narayan (1991) and Paczyhski (1991) investigated the nature of the 
torque for a cold thin disk, but the case of a hot flow has not been studied, (iii) For similar mass 
accretion rates (in Eddington units), the luminosity of a NS accreting via an ADAF is likely to 
be much higher than that of a BH because a NS has a surface while a BH has an event horizon 
(Narayan & Yi 1995b; Narayan, Garcia k. McClintock 1997; Menou et al. 1999). (iv) The spectra 
are expected to be different. 

We discuss in this paper the structure of a hot accretion flow around a NS, or any other 
relativistic star with a surface. The flow under consideration is global and extends radially to very 
large distances (at least thousands of NS radii or more) where it matches onto appropriate outer 
boundary conditions. We do not attempt a detailed analysis of the boundary layer region near the 
NS, where the accretion flow meets the star. We present only an approximate analysis of this region, 
which extends at most a few NS radii above the stellar surface (for a more detailed discussion of the 
physics of the boundary layer see Popham & Narayan 1991; Paczyhski 1991; Titarchuk, Lapidus, 
k Muslimov 1998; Titarchuk k Osherovich 1999). 

The paper is organized as follows. We show in §2 that there is a radially extended region 
around the NS where the flow "settles" with a radial velocity much less than the local free-fall 
velocity. We obtain a self-similar solution for this settling region and show that, surprisingly, the 
density and temperature of this zone are independent of the mass accretion rate. In §3, we discuss 
physical properties of the accretion flow. In §4 we compare the analytical results with numerical 
computations and in §5 we discuss the relationship between the settling flow and an ADAF. We 
conclude with a discussion in S6. 



2. Self-Similar Settling Solution 

We consider a steady, rotating, axisymmetric, quasi-spherical, two-temperature accretion flow 
onto a star with a surface. We use the height-integrated form of the viscous hydrodynamic equations 
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(Ichimaru 1977; Abramowicz et al. 1988; Paczynski 1991; Narayan & Yi 1994): 

M = AttR^pv, (1) 
„^=P^-.l)H-i^(p^), (2) 

where M is the mass accretion rate, R is the spherical radius, p is the mass density of the accreting 
gas, V is the radial infall velocity, U is the angular velocity, ^xiR) = {GM / R^Y^"^ is the Keplerian 
angular velocity, = (?{0p + OerUe/mp) is the square of the isothermal sound speed, Tp^e are the 
temperatures of protons and electrons, 9p^e = kBTp^e/mp^ef^ are the corresponding dimensionless 
temperatures, a is the Shakura-Sunyaev viscosity parameter, J is the rate of accretion of angular 
momentum, Sp and Sg are the specific entropies of the proton and electron fluids, p^ ~ {me/mp)p is 
the mass density of the electron fluid, 7^ and 7e are the adiabatic indices of protons and electrons 
(which, in general, may be functions of Tp and Tg), and q'^,q~, and qcou\ are the viscous heating rate, 
radiative cooling rate, and energy transfer rate from protons to electrons via Coulomb collisions, 
per unit mass. We have assumed that a fraction 6 of the viscous heat goes into electrons and a 
fraction 1 — 5 into protons; it is usually assumed in ADAF models that (5 <C 1, but our analysis is 
general for any value of 5 between and 1. Equations (l)-(5) describe the conservation of mass, 
radial momentum, angular momentum, proton energy and electron energy, respectively. 

For simplicity, we have assumed in equation (1) that the flow is spherical. A more accurate 
treatment would replace R^ with RH, where the scale height H = Cg/Ofc. This would introduce 
minor differences in some of the quantitative results. 

In the case of accretion onto a NS we expect the flow to slow down as it settles on the stellar 
surface, and we expect the density in this settling zone to be significantly higher than for a BH. 
The increased density would cause more efficient transfer of energy from protons to electrons via 
Coulomb collisions and more efficient radiation from the electrons. As we show below, this leads 
to a flow in which g"*", q~ and gcoul are all of the same order, which is very different from the 
case of a BH ADAF, where q'^ ^ q^ , t/coul- Another feature of the settling zone, again the 
result of the large density, is that optically thin bremsstrahlung cooling (which is sensitive to p) 
dominates over self-absorbed synchrotron cooling. We therefore neglect synchrotron emission in 
our analysis. For simplicity, we neglect also thermal conduction. Comptonization is important over 
part of the settling zone. However, we have not been able to derive useful analytical results with 
both bremsstrahlung and Comptonization included. Therefore, we neglect Comptonization in this 
section, and discuss its effects in §3.3. 
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The set of equations (l)-(5) must satisfy certain boundary conditions at the neutron star. 

First, as the flow approaches the surface of the star at i? = RnSj the radial velocity must become 
very much smaller than the local free-fall velocity. Second, the angular velocity must approach the 
angular velocity of the star 0,ns- We use the dimensionless parameter 

^k{Rns) 

to represent the spin of the NS. 

The radius of the star, and its spin, are the two principal boundary conditions applied at 
the inner edge of the accretion flow. We assume that the star is unmagnetized, so there are no 
magnetospheric effects to consider. Two outer boundary conditions, namely the temperature and 
angular velocity of the gas, are determined by the properties of the gas as it is introduced into the 
accretion flow on the outside (e.g. from an ambient medium or from a different type of accretion 
flow such as a thin disk). These outer boundary conditions have little effect on the interior of 
the flow (see §5 of this paper, and Narayan, Kato &: Honma 1997; but see also Yuan 1999). An 
additional important boundary condition is the mass accretion rate M, which is determined by 
external conditions and which we take to be constant. 



2.1. Inner Settling Solution 

We consider first the inner region of the flow, Rns ^ R ^ where Rs = 2GM/(P is the 

Schwarzchild radius. In this region we expect a two-temperature plasma, with Tp > Tg, in which 
the electrons are relativistic and the protons are non-relativistic: Og ^ 1, 9p <^ 1. The viscous 
heating rate of the gas, the energy transfer rate from the protons to the electrons via Coulomb 
collisions, and the cooling rate of the electrons via bremsstrahlung emission are given by 

- "^(^)^ 

qconi = QcouiP^/, Qcoul = 47rrglnA— (8) 

3 

q~ = Qs,RP%, Qff,R = 48a/^e^> (9) 

where is the fine structure constant, Ve is the classical electron radius, In A ~ 20 is the Coulomb 
logarithm, ~ c^Op, and we have neglected logarithmic corrections to the relativistic free- free 
emissivity. The subscript "R" in Qg^R denotes relativistic bremsstrahlung. 

We now make a number of simplifications in equations (l)-(5). First, we neglect the radial 
velocity term vdv/dR in equation (2). Second, we assume that J dominates over MQ,R^ on the 
right-hand-side of equation (3) and we neglect the latter term. Third, we neglect the entropy terms 
in equations (4), (5); that is, we assume that the heating, cooling and energy transfer terms in 
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these equations dominate over the entropy gradient terms. All of these assumptions are justified a 
posteriori below. The equations then read 

M = AttB?pv, (10) 
{^'-^]<)R=\^{pcl), (11) 

''"'-triR'-'' 

(1 - S)q+ = (lo„„l = «- - q*S, (13) 

where q~^, gcouh and q~ are given by equations (7)-(9). 

It is straightforward to show that the simplified equations (10)-(13) have the following self- 
similar solution, 

p = por-'^, 9p = 9por-'^, 9e = 9eor-^/'^ , 

n = ^^or~^/^ V = vor^, (14) 

where r = R/Rs is a dimensionless radius. The normalization coefficients are uniquely related to 
each other as follows 

Opo = ^(1-^')' (15a) 

= felT^J =112^(1^; ^''■'[t^) ' ^^^'^ 
Oo = J^ifos - 7.19 X 10^ m~^s rad/s, (15c) 



Po 



4Qcoul 



a(l - (5)eeos^ - 8.10 X 10-^m-^a(l - 6)9eos'^ g/cm^ (15d) 



j = -6nac'Rlpo9,os = -^^M^a\l-6)yh^l-sY'- (15f) 

2v3 -\/(^couiVfr,R 

Here, ilxo = ^KiRs) and the dimensionless spin parameter is ,s = 0,{Rs) /^ko = ^NS /^KiRNs) = 
^1{R)/Qk{R), as introduced in equation (6). Recall that ,s is a boundary condition of the problem. 

The range of r over which the solution is valid is determined by the twin requirements that 
the protons be non-relativistic and that the electrons be relativistic. The former condition is 
satisfied for any r > 1, while the latter condition requires r < c± 120. A third condition is that 
Comptonization should be negligible (since we have assumed this). As we show in §3.3, this last 
condition requires r > few tens, with the exact limit depending on the NS spin, s, and viscosity, a. 

The radial velocity of the solution is independent of r, whereas the local free-fall velocity 
varies as u// = c/\/2r. Thus, v/vff decreases with decreasing radius. This shows that the solution 
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corresponds to a settling flow and that it is quite different from self-similar ADAFs around BHs, 
where v/vff either is constant (Narayan & Yi 1994, 1995a; Manmoto et al. 2000) or increases with 
decreasing r (Narayan ct al. 2000; Quataert & Gruzinov 2000). Although v/vjf is quite small as 
the flow approaches the NS surface, v itself is still fairly large. At the NS surface, v must reduce 
substantially from its self-similar value. As the numerical results of §4 show, this happens in a 
boundary layer where the accreting material cools catastrophically to a temperature that is orders 
of magnitude below virial. The boundary layer is distinct from the settling zone which is described 
by the above self-similar solution. 

The angular velocity of the gas is a fixed fraction of the local Keplerian angular velocity, the 
ratio being determined by the dimensionless spin s of the star. The gas in the settling solution 
radiates most of the energy dissipated through viscosity. In fact, the rates of viscous heating. 
Coulomb energy transfer and radiative emission are all equal, which is achieved by a suitable choice 
of the density, electron temperature and proton temperature in the gas. The two temperatures have 
universal forms, with only a weak dependence on s, while the density has a strong dependence on 
s. 

The most outstanding feature of the self-similar solution is that, except for the radial velocity, 
none of the other gas parameters has any dependence on M. 

The fact that J is negative implies that the accretion flow removes angular momentum from 
the star and spins it down. This behavior is quite different from that seen in thin disks (Popham & 
Narayan 1991; Paczyhski 1991), where for most choices of the stellar spin parameter s, the accretion 
disk spins up the star with a torque Jthin ^ MflK{RNs)R%s- O^^Y when s is very close to unity 
does the torque become negative. In contrast, for the self-similar solution derived here, the torque 
is negative for all values of s (except extremely small values, see the discussion in §4). Moreover, 
J is independent of M. Equivalently, the dimensionless torque, j = J /M^k{Rns)R%S' which is 
~ 1 under most conditions for a thin disk, here takes on the value 



V^rnp ri _ ^-l/2^ _ ^2^-1^3 _ ^2n3/2 



: r 

„2n3/2 



~ -A3a''m-^s^l-sY' , (16) 

where m = M / M^dd is the mass accretion rate in Eddington units, with MEdd = 1-39 x 10"^*m g/s, 
(for a nominal tj = 0.1), rjsfs ~ 3. Note that —j could be very large at low m. 

We now check under what conditions the approximations we made earlier arc valid. First, we 
neglected the term vdv/dR in equation (2). This is obviously valid since the self-similar solution 
has V = constant. 

Second, we assumed that \J\ » \MVlR^\ in equation (3). Using (16) this condition may be 
cast into the form |j| 2> sr^/^/\/3, or equivalently. 



(17) 
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For a ~ 0.1, s ~ 0.3, and assuming a radial extent of r ~ 10^ for the flow, we require m < 7 x 10~^. 
A direct numerical simulation (§4) shows that the analytical solution is valid even for values of rh 
that are a factor of a few larger than this limit. 

Third, we neglected the entropy terms pp^evTp^edsp^e/ dR in equations (4) and (5). For the pro- 
tons, using the solution (14), we can show that the left-hand-side of equation (4) varies as while 
the right-hand-side varies r~^'^. Thus, with decreasing r, the entropy term becomes progressively 
less important than the other terms, thereby confirming the validity of the approximation. In the 
case of the electrons, the entropy is always small since Sg ~ {meTe/mpTp)sp <^ Sp. 



2.2. Outer Settling Solution 

For r > 10^'^, both protons and electrons are non-relativistic and the solution described in the 
previous subsection is not valid. Interestingly, another self-similar solution may be derived for this 
region of the flow. This solution has a nearly one-temperature plasma with Tp — Tf,<^ Tp, Tg. The 
free-free cooling takes the form 

3 

9~ ~ Qfr,NRp'^e^', Qfr,NR = 5V27r-3/2^/aT^, (9') 

nip 

where cry is the Thompson cross-section, and the subscript NR stands for non-relativistic. Since 
the gas is effectively one-temperature. Equation (13) simplifles to 

q+ ~ q-. (130 

Wc do not need to consider the Coulomb transfer rate gcouh since this quantity is proportional 
to (Tp — Tg) and can be adjusted to have the right magnitude with small adjustments of the two 
temperatures. In the non-relativistic regime, 0e < 1, ~ ("i-e/"ip)6*e 1) and the Coulomb 
transfer rate is 

3 me arc 2^^^ ~ kTe 

^From the condition gcoui — Q~ it follows that 

kTp - kTe ~ i^^^?^c202g-i/^% (19) 

which is exponentially small for 9e < 1- 

The flow is again described by the self-similar solution (14), with the following two exceptions: 

eOr-\ eeO = ^epO, (15b') 



m, 



e 



P = Por , Po = ^77^ ^ = 0.12m OA as gem, 15d') 
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where we have used the fact that the total pressure p = Pp + Pe = "^Pp- 
This solution is valid only if m is quite small, cf. equation (17): 

m < 2.2 X l^~^Oil isl 2,rl^''^, (20) 

where rs = r/10^, ao.3 = a/0.3, and sq.i = s/0.1. For greater m, a self-similar power-law solution 
does not exist; numerically computed solutions exhibit non-power-law behavior, as discussed in §4. 



3. Properties of the self-similar solution 
3.1. Spin- Up/Spin-Down of the Neutron Star 

The rate of spin-up of the accreting NS is given by 

j^{Ins^) = j-Mn{RNs)R%s 

~ -A3s^a^MEdd^K{RNs)R%S: (21) 

where Ijys is the moment of inertia of the NS. We have made use of the fact that > \Mn{RNs)R%s\ 
for the self-similar solution, and used equation (16) for J. The negative sign in the final expression 
implies that the accretion flow spins down the star. The above equation is for an unmagnetized NS. 
If the NS has a magnetosphere, the inner edge of the accretion flow is at the magnetospheric radius, 
Rm- In this case, let us define s by ^ns = sUK{Rm) = s^K{RNs){Rm/ Rns)"^^"^ ■ Substituting 
this in equation (21) with Ins = constant and integrating, we obtain 



50 _ Ins f Rm \~^'^ 



VTTiT^' 86s2a2MEdd-R^5 \RnsJ 



(22) 



where sq = s(t = 0). The same result is valid for an unmagnetized NS by setting R^ = Rns- The 
quantity r is the characteristic spin-down time of the NS. For a spherical NS of constant density. 
Ins = 2MR'^g/5 = (0.8 x 10^^ g)mR%g. Substituting this expression, we obtain the spin-down 
rate Pns/Pns = with 

r 6.7 X lO^^.'^a'^ (|^) s = 2 x 10«.oi«oi (|^) yr. (23) 

Note the remarkable fact that the spin-down time scale is independent of the mass of the NS, and 
the mass accretion rate! For the magnetic case, the rate depends on the radius ratio Rm/RNS- 

It is customary to express the spin-down rate as Pns/Pns- Writing 
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and ^k{Rns) — lO^m^ \ rad/s, where Rns = ^Rs and mi.4 = M/{1AMq), we obtain 

^ ~ 2.4 X IQ-^V^^a^s^ s"^ = 2.7 x IQ-^^m^^a^.gS^.s s"^ (25) 

where so.5 = s/0.5. This spin-down rate is in good agreement with observational data on the spin- 
down of X-ray pulsars for which Yi, Wheeler Sz Vishniac (1997) invoked ADAFs: 4U 1626-67 has 
P/p2 8 X 10-^3 s-2 and P = 7.7 s; OAO 1657-415 has P/P^ ^2x 10" g-2 p = 38 s, and 
GX 1+4 has P/p2 « 3.7 X 10-^2 s-2 p ^ s. Since the spin-down rate is quite sensitive 
to a and s, the observed data in individual systems can be fitted by small adjustment of these 
parameters. 



3.2. Luminosity 

In computing the luminosity of the accretion flow, we must allow for the energy release in both 
the boundary layer and the self-similar settling zone. We calculate their luminosities separately. 

Radiation from the self-similar settling flow may be calculated following the methods described 
by Popham &; Narayan (1995) for a thin disk. This method assumes that the luminosity at a given 
radius is determined by the local viscous energy production. This is a legitimate approximation for 
the settling flow in which q" = q'^. Keeping only the dominant terms, we find 

GMnsM / 1 2 \ ,v Z"^""* dP 
Lss= {l + -s^-3sj+M j^^ — , (26) 



R-in V 2 / JPin P 

where Rm = Rns + ^bl is the inner radius of the self-similar zone and Abl Rns is the thickness 
of the boundary layer. Here the first two terms, (1 + s^/2), represent the luminosity, associated 
with potential energy of the infalling gas, Lpot, the third term —js is the luminosity, associated 
with the rotational energy extracted from the star, Lj-ot (note, j < in the self-similar solution), 
and the final integral is the "enthalpy correction", -Lenth- Using the analytical solution (14)-(15) 
and assuming Pout = for simplicity, we obtain 

2 rn f ^ 



Lpot = MEddc'^^ (^^ + 2" ) ' ^^^^^ 



a2 



Lrot = 43MEddc';^ il-S)^/^s^{l-s^), (27b) 

2rNS 

Lenth = -MEddC^T^{l-S^). (27c) 

2rNS 

Note that the leading terms in Lpot and Lenth cancel each other exactly. The luminosity of the 
self-similar settling zone is thus 

Lss ^ 6.2 X lO^'^mrg ^m_2So.i + 8.9 x lO^^mrg ^a^^s^ ^ ergs-\ (28) 



-10- 



where m_2 = m/0.01, sq.i = s/0.1, rs = r^rs/S, and we have assumed s <C 1. Note that luminosity 

not associated with with dissipation of rotational energy, represented by the first term in equation 
(28), is much less than the commonly assumed GMjsfsM /Rj^s- This is because the negative 
enthalpy term has large magnitude, as a result of the fact that the settling flow is akin to a 
pressure supported, quasi-stationary atmosphere. 

The second term in equation (28) is the luminosity of the settling zone. Since the self-similar 
solution for this zone is independent of rh, the luminosity too shows no m dependence. Indeed, 
the luminosity remains finite even as m — 0. How is this possible, and where does the energy 
come from? The answer is that the luminosity of the settling zone is supplied by the central star. 
As the star spins down, it does work on the accretion flow and the energy released comes out as 
bremsstrahlung radiation. 

The boundary layer luminosity requires a different method of calculation since viscous energy 
production is negligible in this zone: Q ^constant, and so g+ oc {dVL/dR)^ ~ 0. As the accreting 
gas cools in the boundary layer, starting from a nearly virial temperature ~ 10^^ K on the outside 
down to the NS temperature ~ 10^ K near the surface, the thermal energy in the gas is emitted 
as radiation. To estimate the luminosity, we use the energy balance equation, which is the sum of 
equations (4), (5): 

pv dc^ r, dp 7 dc^- dP 

-q~ = — ^--vcl— = — —pv—-v . (29) 

^ -f-ldR "dR -i-V dR. dR ^ ' 

We can neglect the dP/dR term because the pressure P is essentially constant in the boundary 
layer. To obtain the luminosity we integrate over the boundary layer 

Lbl = ! q~4TrR^dR = - [ -^^AnR^ pv^ dR = -^^MAct. (30) 
J J7-1 dR 7-1 

Since cj starts from nearly virial value and reaches close to zero, Ac^ ~ GMns/Rns- More 
precisely, Ac^ = c^A{9p + 6^) — c^Op^rJ^^g. Therefore, the boundary layer luminosity is 

Lbl = ^^MEddC^^^ (l - s^) ~ 1-7 x lQ^^mm-2r^^ , (31) 
7-1 QrNS 

where we have assumed 7 = 5/3. The total luminosity of the system is L = L55 -|- Lbl- 
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3.3. Effect of Comptonization 

Using the self-similar solution (14), (15), we may readily estimate the electron scattering optical 
depth and the y-parameter.^ The optical depth is 

~ ao.isl,ir~^, (32) 

where Kes = OT/fnp is the electron scattering opacity for ionized hydrogen. Since r > 3, we see 
that Tes < 1/3 for reasonable parameters and the radiation is optically thin to electron scattering. 
The y-parameter is 

y = 16^,Ves-2xl0Vl-'J)"''''(l-^')'/Vr-2 

~ 2 X lO^ao.iSo.i^"^- (33) 

The radius at which y 1 is 

re~45a;/i^so.i. (34) 

Above this radius the inverse Compton scattering is small and the self-similar solution is valid. For 
r < Tc, however, Comptonization is important and the electron temperature profile will be modified 
from the self-similar form. Since the electron-proton coUisions are relatively weak (the plasma is 
two-temperature), other quantities, e.g., the density, proton temperature, etc., are unaffected. 
Comptonization is unimportant for low-viscosity flows, a < 0.01 around slowly rotating NSs, s < 
0.01, because then Tc < r^s- 



3.4. Spectrum 

We now estimate the spectrum of radiation emitted from the settling accretion flow. Let us 
neglect inverse Compton scattering for the moment. The relativistic bremsstrahlung emissivity 
is approximated as Ci, oc exp — {hv/kTi,) erg cm~'^ s^^ Hz~^. Therefore the luminosity per unit 
frequency is 

/•oo 

oc / p^e-^'''^'^''2'KB?dR 

J Rns 
/•oo 

OC / r^e-''^dt(xi^'^T{-2,iy/iym), (35) 

Jl/Um 



^Here we just estimate where the effect of Comptonization becomes significant. For better analytical approxi- 
mations see, for instance, Dermer, Liang, & Canfield (1991); Titarchuk & Lyubarskij (1995). (In the latter paper, 
the expression for y is not given, but it can be inferred using equation [24]: y = Tes [(a -|- 3)0/(1 + 0) + AdQ^^O^].) 
Comptonization of free-free radiation has also been considered by Titarchuk (1989). 
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where T{a,z) = t"'~^e~*dt is the incomplete gamma-function and u^n = kTe{RNs)/h is the 
maximum frequency. Above Um the spectrum falls exponentially and below Um it is nearly flat. We 
may, thus, replace the exponential in the integral with a square function which is equal to unity 
for v < and for ly > u^. With this approximation 

3 Lss (, v^- \ „„, 

where Lss = / Lydu is the total luminosity of the self-similar flow, represented by equation (28). 
The break frequency, f^, is roughly given by hvrn ~ 2.7 MeV for a typical electron temperature 
7'e,max ~ lO''^'^-^ °K [cf., equation (15)]. At a typical x-ray energy, /iz/ ~ 3 keV, the observed 
luminosity per decade is 

z^L, ~ 1.7 X lO^ima^ (^^^^ erg s'^, (37) 

i.e., vLi, ~ 1.5 X 10^^ for a 300 Hz neutron star (sq.i ~ 1-6). The luminosity per decade is much 
greater at higher photon energies and may be as high as ~ few x 10^^ — 10^^ erg/s at hv ~ MeV. 

As shown in the previous section, Comptonization becomes important below the radius Vc- At 
Tc, y ~ 1 and the electron temperature is 

re(rc) ~ 2.7 MeV/^ ~ AQQa^]'^ s^^^"^ kcV. (38) 

For r < Tc, the electron temperature will be determined self-consistently by Compton cooling rather 
than by bremsstrahlung emission. Computing the spectrum from this region is beyond the scope 
of the paper. We also do not attempt to calculate the spectrum of the radiation from the boundary 
layer. 



3.5. Bernoulli parameter 

It is known that the Bernoulli parameter of the accreting gas in BH ADAFs is positive for a wide 
range of r (Narayan k, Yi 1994, 1995a; Narayan, Kato k, Honma 1997), and it has been suggested 
that the positive Bernoulli parameter may trigger strong winds or jets in these systems (Narayan 
& Yi 1994, 1995a; Blandford & Bcgclman 1999, but see Abramowicz et al. 2000). Igumenshchev &: 
Abramowicz (1999, 2000) confirmed with numerical simulations that strong outflows are produced 
from BH ADAFs when a ~ 1. 

Normalizing the Bernoulli parameter, Be, by (fix-R)^) and using equations (14), (15), we find 
that the self-similar settling flow has 

c2^2 ^7-1 3 
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(39) 



3(7-1) 2 3(7-1)' 

where 7 is the mean adiabatic index of the flow and in the last expression we have neglected the 
term in Vq since it is negligible deep inside the settling region. 

We see that the second term in the final expression is always negative. This term is proportional 
to s^, which means that a more rapidly spinning NS stabilizes the accretion flow against outflows 
more effectively than a slower spinning star. The physical explanation is as follows. The centrifugal 
force increases with increasing s, which has the effect of making it easier for the gas to escape. At 
the same time, however, the centrifugal force causes the radial infall velocity to decrease, which 
increases the time available for cooling. The temperature and the gas pressure go down, making the 
gas more gravitationally bound. The net contribution of the two effects turns out to be negative. 

The first term in the last line of equation (39) can be either positive or negative, depending on 
the value of 7. Combining the two terms, we find that the gas is gravitationally bound and unable 
to flow out in a wind (i.e. 6 < 0) if the adiabatic index satisfies 



7 > -bl^- ^4°) 



For <C 1, the condition is 7 > 1.5; that is, the the accretion flow can produce a wind and/or a 
collimated outflow only if 7 < 1.5 and is stable to such outflows if 7 > 1.5. Normally, we expect 7 
to be close to 5/3 for the accreting gas. 



3.6. Stability to Convection 

It is well known that if the entropy increases inwards in a gravitationally-bound non-rotating 

system, the gas is convectively unstable; otherwise the flow is stable. The specific entropy profile 
in the settling accretion flow around a NS can be readily calculated from equations (4), (5) using 
(14), (15). This gives 

ds _ k 1 d / c2 \ _ 27 - 3 1 

We see that the entropy increases outwards for 7 > 1.5 and inwards for 7 < 1.5. Hence if 7 > 1.5 
the flow is stable against convection, while if 7 < 1.5 the flow is convectively unstable. 

In the presence of rotation, the analysis is a little more complicated. Narayan et al. (2000) and 
Quataert & Gruzinov (2000) discuss the generalization of the Schwarzchild criterion for accretion 
flows with rotation. If the gas motions are restricted to the equatorial plane of a height-integrated 
flow, convective stability requires the following effective frequency to be positive: 



(42) 
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where N is the Brunt- Vaisala frequency and k is the epicycHc frequency, «; = O for oc For 
a power-law flow with p oc and Q{R) = sCl^ ^ R^^^"^ with = 1 — (1 + a)cQ, this criterion 
may be written as follows (see Narayan et al. 2000 for more discussion) 




Since for the self-similar settling solution a = 2, the stability criterion (43) becomes 




[(27-3) + s2(3-7)] >0 



(44) 



which yields that the flow is convectively stable if 



This condition is different from the stability criterion against outflows, given in equation (40). 

Following the techniques developed by Quataert & Gruzinov (2000), Narayan et al. (2000) 
have also presented a more general analysis of convection in a self-similar accretion flow. This 
analysis, which docs not restrict motions to lie in the equatorial plane, assumes that t;^ and 
are independent of the polar angle 9 (as is valid for a marginally convectively stable system, cf 
Quataert &; Gruzinov 2000) . Narayan et al. (2000) find that the most unstable region of the flow is 
near the rotation axis, = 0, tt. They show that the marginal stability criterion for this polar fluid 
coincides with the condition for the positivity of the Bernoulli parameter. That is, a flow which is 
convectively stable at all 6 has a negative Bernoulli parameter, while a flow which is convectively 
unstable for at least some values of has a positive Bernoulli parameter. (The Bernoulli parameter 
itself is independent ol 9.) 

We have verified this result for the solutions presented in this paper. Specifically, when we 
apply to our solution the more general convective stability criterion given by equation (A9) of 
Narayan et al. (2000), we recover the condition (40) above, namely that the self-similar flow is 
convectively stable if and only if 



We have numerically solved the system of height-integrated two-temperature fluid equations 
(l)-(5) with boundary conditions. In the energy equations we assume that viscous dissipation only 
heats the protons (i.e. (5 = 0). We include energy transfer from protons to electrons via Coulomb 
collisions, and we take the cooling of electrons to be purely by free-free emission, as discussed in 
§2. For these processes we use the expressions given in Narayan &; Yi (1995b), which smoothly 




(46) 



4. Comparison with Numerical Results 
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interpolate between the regimes of non-relativistic and relativistic electrons. We take into account 

the variation of the electron adiabatic index 7e with temperature (Chandrasekhar 1939; Esin et al. 
1997), using a simple interpolation formula from Gammie &; Popham (1998). We assume that the 
protons have 7p = 5/3. 

We employ the gravitational potential of Paczyhski & Wiita (1980) to mimic the effect of 
strong gravity near the NS surface. In this potential the Keplerian angular velocity takes the form 

Note that the analytical work presented in the previous sections is based on a Newtonian potential. 

We specify the boundary conditions as follows. We take the outer boundary of the flow to 
be at Tout = 10^. At this radius we specify that the angular velocity is equal to its value in the 
self-similar ADAF solution of Narayan & Yi (1994), and that the proton and electron temperatures 
are both equal to the self-similar ADAF temperature. We assume that the accreting star is a IMq 
neutron star with a radius Rns = 3i?s = 8.85 km, unless stated otherwise. At i? = Rns, we specify 
the value of the NS spin parameter, s = ^ns/^k{Rns), and we require the proton temperature 
of the flow to be T = few x 10^ K ^ Tyij-i^i- (We do not assume that the electron and proton 
temperatures are equal, but in fact they are equal.) We do not constrain the density of the gas in 
any way at either boundary. 

The numerical problem as posed here has a family of solutions characterized by three dimen- 
sionless parameters: the mass accretion rate m (in Eddington units), the NS spin s (in units of 
the Keplerian angular velocity at the NS surface), and the viscosity parameter a. The angular 
momentum flux J, or equivalently the dimensionless flux j = J/MQ,k{Rns)R'ns' ^ eigenvalue 
of the problem. 

Figure 1 shows representative solutions for a = 0.1 and a range of values of m and s. The 
solutions clearly have three radial zones. For r > 10^-^, there is a one-temperature zone in which the 
gas properties vary roughly as power-laws of the radius. For r < 10^-^, there is a second power-law 
zone with a two-temperature structure. Finally, close to the NS, the flow has a boundary layer 
region. In this final region, the gas experiences run-away cooling, the velocity falls precipitously, 
and the density increases very rapidly. This region of the flow does not have power-law behavior. 

The numerical solutions are unreliable in the boundary layer; thermal conduction (not included 
in the calculations) is probably very important here, and optical depth effects (also not included) 
will modify the radiation properties significantly. The solutions are suspect also in the inner region 
of the two-temperature power-law zone, where Comptonization is likely to be important. Outside 
these regions, however, the numerical solution is expected to be accurate. 

The numerical results in the two-temperature power-law zone below r ~ 10^^^ agree quite well 
with the analytical solution presented in §2.1. Curves corresponding to a given value of s and 
different values of m coincide with one other to very good accuracy, as predicted by the analytical 
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solution. This is best seen in the profiles of p and Changing s causes an up/down shift of the 

curves but docs not affect the slopes of the curves. The temperature profiles are sensitive to the 
spin s, especially for large values of s. The radial velocity varies approximately as v oc m and is 
roughly consistent with u oc for s > 0.1. 

For r > 10^'^, the numerical solutions are in reasonable agreement with the one-temperature 
self-similar solution described in §2.2. The agreement is less perfect than in the previous zone. This 
is primarily because the analytical solution requires a very low value of m in order to be valid as 
far out as the outer radius r ~ 10^ [cf., equation (20)]. The numerical models shown have larger 
values of m than this limit. 

As we discussed in §2 and §3, the transport of angular momentum in a hot settling flow differs 
dramatically from the well-known behavior of a thin disk. The solid line in Figure 2 indicates the 
dependence of the dimensionless angular momentum eigenvalue j as a function of the dimensionless 
NS spin s. The long-dashed line indicates the corresponding results for a thin disk (Popham & 
Narayan 1991), where j ~ +1 for most values of s, and goes negative only for stars nearly at 
break-up. In contrast, in the settling flow, j is negative for almost all values of s. For the particular 
choice of parameters, namely a = 0.1, rh = 0.01, s ~ few x 10~^, we find that j ^ — few. 

The short-dashed line in Figure 2 shows the analytical formula for j, as given in equation 
(16). The agreement with the numerical results is good for a wide range of s below about 0.5. For 
s > 0.5, the numerically determined j levels off at a constant negative value, whereas the analytical 
result shows \j\ decreasing rapidly. The main reason for the discrepancy is the neglect of the ram 
pressure term in the radial momentum equation in the analytical work. Note the interesting fact 
that super-Keplerian accretion (s > 1) is, in principle, possible (provided one can arrange to have a 
star with super-Keplerian rotation) . In a super-Keplerian flow, the ram pressure of the infalling gas 
supplies the radial momentum needed to push the gas onto the NS. For extremely small s <^ 0.1, 
the numerical solutions show j to be slightly positive; we find j ~ 10~^, as indicated in the lower 
panel in Figure 2. Here again the self-similar solution, which predicts a small negative value for 
j, breaks down. For very small s, the J term in equation (3) is comparable to or smaller than 
the MQ.B? term which was omitted in deriving the analytical solution. This is the reason for 
the discrepancy. The precise value of s at which the analytical solution breaks down depends on 
the choice of parameters (a, m), but is essentially independent of the outer radius. Thus, the 
transport of angular momentum through the flow and the spin-down of the star axe determined by 
the boundary conditions at the stellar surface, and not by the outer boundary of the flow. 

One of the most surprising features of the self-similar solution is that the angular momentum 

flux J is independent of m; equivalently, the dimensionless eigenvalue j is oc m~^. The solid line in 
Figure 3 is a plot of — j as a function of m for a flow with a = 0.1 and ,s = 0.3, as determined from 
the numerical solutions. The dashed line indicates for comparison the analytical scaling, j oc m~^, 
with the proportionality constant given in equation (16). Note the very good agreement between 
the numerical results and the analytical self-similar solution. 
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The highest value of m up to which we could obtain a numerical solution is merit = 0.0313. 
Beyond this critical value, there is no hot solution. (The value of rhcrit depends on s, a and rout-) 
For larger m, the density in the flow is so high that there is runaway free- free cooling and the gas 
is unable to remain hot. We presume that the accretion then occurs via a thin accretion disk. 

Figure 4 plots, for selected values of m, the luminosity per logarithmic interval of D, where D 
is the fractional distance from the NS surface: 

i5=<^^. (48) 

Rns 

We see that the luminosity at the peak of the curve is very insensitive to m. The emission in the 
peak corresponds to radiation from the settling flow. This emission represents energy released by 
the spin-down of the NS, and its luminosity is independent of rh [cf. equation (28)]. For radii below 
the peak, the curves do show a dependence on rh. The radiation here corresponds to boundary 
layer emission, which is proportional to rh according to equation (31). 

All the models described above have Rns = ^Rs- However, different equations of state predict 
slightly different NS radii, Rns = 2 — 4Rs. We have computed numerical models for this range of 

Rns and we find that the subsonic settling solution exists for the whole range. Qualitatively, the 
solutions with different Rns are very similar. As Rns decreases, the peak temperature is higher, 
as expected for the deeper potential. 



5. Relationship of the Settling Flow to an ADAF 

The accretion solution we have discussed so far radiates all the energy dissipated by viscosity, 
and is therefore "cooling-dominated." On the other hand, it is known that a hot flow around a black 
hole is an "advection-dominated" accretion flow (ADAF). Both solutions have a two-temperature 
structure for r < 10^'^ and both are very hot (nearly virial) for all r. How are these two types of 
accretion flows related to each other? 

By solving the equations numerically for different boundary conditions, we have found that 
the two solutions are part of a single sequence of solutions in which the spin of the star, s, plays a 
pivotal role as a control parameter. For relatively rapidly rotating stars, with s > 0.1, we obtain 
the settling solution in our numerical experiments. However, as s is decreased, we find that the 
settling solution smoothly transforms to an ADAF-type solution, which becomes well-established 
for s < 0.01. The transition is not sharp, so it is difficult to identify a specific transition point 
s = St at which the transformation occurs. Numerical experiments indicate that the value of st 
(however it is defined) is not very sensitive to Rout-,1-, and rh and is, roughly, st ~ 0.04 — 0.06. 

The change of the nature of the fiow as s is varied is illustrated in Fig. 5. The solid and dotted 

curves correspond to two solutions with s = 0.3 and s = 0.01, respectively, with all other boundary 
conditions being the same. We see that the solutions are markedly different from each other. This 
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is most clearly seen in the profiles of density, where the s = 0.3 model has a logarithmic slope 

of -2, as appropriate for the cooling-dominated settling solution described in this paper, and the 
s = 0.01 model has a slope of -3/2, as expected for a standard self-similar ADAF (Narayan & Yi 
1994, 1995a). There is a similar difference also in the profiles of the radial velocity, where the two 
solutions have logarithmic slopes of -1/2 and 0, respectively. 

An interesting feature of the s = 0.01 ADAF-type solution is that it consists of two distinct 
segments. For large radii (in Fig. 5, for radii outside r ~ 20), the flow corresponds to the standard 
ADAF discussed in the literature, with the scalings 

p oc r"^/^ cl oc r"S n oc r"^/^ v oc r"^/^. (49) 

However, at smaller radii, the numerical solution indicates the presence of a second advection- 

dominated zone, a "settling ADAF," which was first seen in numerical calculations described in 
Narayan & Yi (1994). This settling ADAF is seen in Fig. 5 as a zone that lies between the boundary 
layer region and the outer standard ADAF, with different slopes for p and v. The radial extent of 
the settling ADAF zone may be quite large and, in general, depends on the flow parameters and 
boundary conditions. 

A self-similar model of the settling ADAF may be readily obtained as follows. In an ADAF, 
energy is not radiated, therefore q~ = 0. Close to the star U ~ constant, therefore g"*" = 0. 
Equations (4), (5) then simplify to the condition of entropy conservation, ds/dR = 0, which yields 
oc p^~^. As the material settles on the star, its radial velocity decreases and we have v <^ 
Vff, <C ^K- Then, from equation (2), it follows that the temperature of the gas is nearly virial. 
Other quantities are determined straightforwardly, so that we have 

1 2 -1 27-3 

p OC rf-1 , oc r , ~ const., v oc r t-i . (50) 

The infall velocity decreases with radius if 7 < 1.5, and increases if 7 > 1.5. To highlight the 
difference between the standard ADAF and the settling ADAF, we have chosen 7 = 4/3 in the 
solutions shown in Fig. 5. 

Finally, the long-dashed curves in Fig. 5 correspond to a solution with s = 0.3 for which we 
have increased the outer boundary value of T by a factor of 10. We see that the solution in the 
interior is not sensitive to the outer boundary conditions (within a reasonable range, of course). 

6. Summary and Discussion 

In this paper, we have presented analytical and numerical solutions that describe a hot, viscous, 
two-temperature accretion flow onto a neutron star (NS). The results are relevant also for accretion 
onto other compact stars with a surface, e.g. white dwarfs. To our knowledge, this is the first study 
of viscous fluid dynamics for a hot flow around a NS. 
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The presence of a surface modifies the nature of the flow relative to the case of a black hole. 
We show that the accretion flow has an extended settling region in which the radial velocity v is 
constant; v is also small relative to the local free-fall velocity. The density in the settling region 
varies as p oc r~^, and the angular velocity has a Keplerian scaling, O = s^Ik oc r~^/^, with s being 
a constant. Here, r is the radius in Schwarzchild units, and the value of s is set by the spin of the 
NS: s = ^Ns/^K{rNs)- At the inner edge of the settling region, there is a narrow boundary layer 
in which the velocity falls extremely rapidly and the density increases sharply to match the surface 
density of the NS. 

The settling region consists of two distinct zones. In the inner zone, r < 10^-^, the gas is 
two-temperature, with the proton temperature varying as Tp oc and the electron temperature 
varying as Tg oc r~^/^. These scalings are derived assuming that electrons radiate primarily by 
free-free emission and that energy transfer from protons to electrons occurs via Coulomb collisions. 
We have derived a completely general analytical self-similar solution for this region which agrees 
well with numerical results. 

In the outer zone of the settling region, r > 10^-^, the gas is one-temperature, Tp Tg oc r~^\ 
here again, we derive an analytical self-similar solution which agrees reasonably well with numerical 
results. 

The most surprising feature of the settling region is that nearly all the gas properties are 
independent of the mass accretion rate m; only the radial velocity shows a dependence: v oc m. 
Since the density and temperature are independent of m, the luminosity is also independent of m. 
Indeed, the settling solution is valid — with a finite luminosity — even in the limit when m ^ 0. 
Clearly, the luminosity does not originate from the gravitational release of energy as mass accretes 
onto the NS. 

The magnitude of the luminosity is very sensitive to the spin parameter s of the NS, varying as 
the fourth power of this quantity, equation (28). For s ~ 0.1, as appropriate for the millisecond X- 
ray pulsar, SAX J1808.4-3658, the model predicts an X-ray (few keV) luminosity vLn ~ 10^^ ergs""^ 
(§3.2). This estimate does not include the contribution from Comptonization in the inner regions 
of the flow, which might increase the X-ray luminosity by an undetermined amount. We note that 
quiescent X-ray luminosities of soft X-ray transients, including SAX J1808. 4-3658, are generally in 
the range vLu < 10^^ ergs~^ (Narayan, Garcia & McClintock 1997; Asai et al. 1999; Menou et al. 
1999; Stella et al. 2000). 

The angular momentum flux J in the settling solution is dominated by the viscous transport 
term rather than the advection term MQ^R^- Consequently, J is negative, i.e. the angular 
momentum flux is oriented outward, and the accretion flow spins down the star. The analytical 
solution predicts that spin-down occurs for all values of the spin parameter s of the central star. 
The numerical solutions by and large confirm this; for plausible parameters, spin-up is seen only 
for extremely small values of the spin parameter, s < 0.005 (cf Fig. 2). (The exact value of s at 
which J changes sign depends on a and m, but is relatively independent of the position of the outer 
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edge of the flow.) This behavior is very different from the case of a thin accretion disk (Popham &: 
Narayan 1991; Paczynski 1991), where one finds that the star is spun-up for nearly all values of s, 
and spin-down occurs only for s close to unity (break-up hmit). 

Another surprising feature of the settling solution is that J, like nearly all other quantities, is 
independent of m. Indeed, the settling zone behaves like a stationary zone (since v is very small), 
and essentially acts like a conventional "brake," slowing down the star by viscosity. The brake can 
operates even ifm^O, so long as there is a static atmosphere of the self-similar form and there is 
a sink for the angular momentum, say an external medium, at large r. Furthermore, the luminosity 
of the settling flow is almost entirely from the energy released by the viscous braking action. That 
is, the luminosity ultimately is fed by the loss of rotational kinetic energy of the star, and not by 
gravity. 

This result has an interesting consequence. In accretion flows around black holes, one defines 
an efficiency factor rj by comparing the accretion luminosity Lace to the rest mass energy of the 
accreting gas, r/ = Lacc/Mc^. It is well-known that r] = 0.06 for a thin accretion disk around a 
Schwarzchild black hole and that the value increases tor) = 0.42 for a maximally rotating Kerr hole. 
A number of interesting ideas have been discussed in the literature for increasing the efficiency of 
an accretion flow around a black hole; these involve tapping the rotational energy of the black hole 
using magnetic fields or viscosity (Blandford & Znajek 1977; Krolik 1999; Gammie 1999). The 
general rclativistic dragging of inertial frames by the spinning hole plays an important role in the 
mechanism of (Blandford & Znajek 1977). 

For the hot settling solution described in this paper, the luminosity is almost entirely from 
the rotational energy of the star. Since Lace is independent of m, the efficiency scales as 77 oc 
and 77 — 00 as m — 0. Thus, it would appear that accretion flows can tap the rotation energy of 
a star with a surface more easily than the energy of a spinning black hole. SibgatuUin &; Sunyaev 
(2000) showed that the extraction of rotational energy of a NS may result in very high boundary 
layer efficiencies, up to ~ 0.67, in the counter-rotating NS-thin disk systems, as well. As in the 
case of the black hole, the energy extraction works best when the star is spinning rapidly: rj cc s^. 
Interestingly, the energy extraction is not a general relativistic effect — our analytical solution is 
based entirely on Newtonian physics. 

Yi, Wheeler & Vishniac (1997) and Yi & Wheeler (1998) recently suggested that the sudden 
torque-reversal events seen in some accretion-powered pulsars may be due to the accretion flow 
switching between a Keplerian thin disk and a hot a sub-Keplerian state akin to an ADAF. Our 
work lends support to this suggestion. We flnd that the torque does reverse in sign between a thin 
accretion disk and a hot settling flow for almost any reasonable stellar spin parameter. We also 
find that the magnitude of the spin-down torque exerted by the settling flow is comparable to the 
measured value of the torque in torque-reversing pulsars (§3.1). 

It is worth emphasizing that the settling solution described here is quite distinct from the 
self-similar ADAF solutions derived for black hole accretion (Narayan &; Yi 1994, 1995a; Honma 
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1996; Kato & Nakamura 1998; Blandford & Begelman 1999; Manmoto et al. 2000). All the black 

hole solutions described in the literature have density varying relatively mildly with radius: p oc 
r~"^/^ — r~^/^. Our settling solution has p oc r~^. Also, the black hole solutions are advection- 
dominated, whereas our solution radiates essentially all the energy it generates through viscous 
dissipation. There are also, as we now discuss, significant differences in the sign of the Bernoulli 
parameter and in the stability to convection. 

Narayan k. Yi (1994, 1995a) showed that their self-similar ADAF solution has a positive 
Bernoulli parameter so long as the adiabatic index 7 of the gas is less than 5/3. They argued on 
the basis of this that ADAFs are likely to have strong outflows and winds (but see Abramowicz 
et al. 2000). Such strong outflows were confirmed with numerical simulations by Igumenshchev & 
Abramowicz (1999, 2000); they found outflows for large values of the viscosity parameter: a ~ 1. 
Blandford &: Begelman (1999) developed a self-similar model with inflow and outflow (the ADIOS 
model). For the settling solution described in this paper, the Bernoulli parameter is positive only 
if 7 is less than 7cr?t, where ^crit = 1-5 for a slowly-spinning star and is smaller than 1.5 for a 
rapidly-spinning star [cf equation (40)]. Since the hot ionized two-temperature gas in the flow 
is likely to have 7 close to 5/3 at most radii, we expect the Bernoulli parameter to be generally 
negative. Therefore, we do not expect a strong outflow. Of course, this conclusion assumes that 
we do not have dynamically important magnetic fields in the flow. 

Narayan &: Yi (1994, 1995a) also showed that their ADAFs are convectively unstable for a wide 
range of conditions. The convective instability has been seen in numerical simulations in which 
the viscosity parameter is assigned a low value: a < 0.1 (Igumenshchev, Chen, &; Abramowicz 
1996; Igumenshchev &; Abramowicz 1999, 2000; Stone, Pringle, & Begelman 1999; Igumenshchev, 
Abramowicz, & Narayan 2000). Self-similar solutions for convection-dominated accretion flows 
(CDAFs) have been derived by Narayan et al. (2000) and Quataert & Gruzinov (2000). For the 
settling solution described in this paper, we find that the gas is convectively unstable only for the 
same low vahics of 7 for which the Bernoulli parameter is positive. (Narayan et al. 2000 showed 
that, in general, for self-similar flows the criterion for the Bernoulli parameter to be positive is the 
same as the criterion for the flow to be convectively unstable, cf. §3.6.) Thus, we do not expect 
hot settling flows around NSs to be convectively unstable, or to have a distinct CDAF mode of 
settling. 

From numerical experiments we have discovered the interesting property that the settling flow 
can continue well inside the last stable orbit (down to at least rjvs ~ 2) and yet remain subsonic at 
all radii. This is despite the fact that the numerical models employ a pseudo-Newtonian potential 
which mimic the last stable orbit for test particles. We find that the structure of the flow is 
qualitatively similar for flows with r^s < and r^s > 3. This is very different from black hole 
ADAFs which become supersonic close to the central object. 

We have also shown (§5) that the settling solution and the ADAF are not two physically 
distinct solutions, but are related to each other. As the neutron star spin is decreased, we flnd that 
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the settling flow smoothly transforms to an ADAF-type solution. The transformation proceeds 
over the spin range s ~ 0.01 — 0.1. 

The settling solution is hot but cooling-dominated. It is thus most closely related to the 
two-temperature solution discovered by Shapiro, Lightman & Eardley (1976). The SLE solution 
is known to be thermally unstable (Piran 1978; Wandel & Liang 199?; Narayan & Yi 1995b), so 
one may wonder about the thermal stability of our settling solution. We defer discussion of this 
important topic to a future paper. 

Could there be non-settling solutions around NSs, and could such flows be more analogous to 
the black hole ADAF, ADIOS, and CDAF solutions? Indeed this is possible if the NS radius is 
small enough. We could imagine, for instance, a standard black hole-like flow around a NS, which 
makes a sonic transition to a supersonic state, and then stops suddenly at a standing shock at 
the surface of the NS. Such a solution would be dynamically consistent, and, except for the shock, 
would be very similar to a black hole flow. However, for such a solution to exist, one requires the 
radius of the NS to be smaller than the sonic radius of the flow. The latter radius is estimated 
to be in the range Vsonic ~ 2 — 5, depending on the value of a (Narayan, Kato & Honma 1997; 
Popham Sz Gammie 1998), whereas NS radii are in the range tats ~ 2 — 4 for typical NS equations 
of state. Thus, for some choices of a and some equations of state, rjv5 could be smaller than 
'''sonic- In such cases, we could have four different hot solutions around a NS: (i) we could have 
the self-similar settling solution of the kind presented in this paper (as we discussed in §4, we find 
subsonic settling solutions for any choice of Rns in the range 2 — 4:Rs)', or (ii) we could have a 
Narayan, Kato & Honma (1997)-like ADAF solution with a shock at the NS surface; or (iii) we 
could have a Blandford &: Begelman (1999)-like ADIOS solution with a shock at the NS surface; 
or (iv) we could have a Narayan et al. (2000) and Quataert & Gruzinov (2000)-like CDAF, again 
with a shock at the NS surface. In the latter three cases, we expect the NS to be spun-up rather 
than spun-down by the accretion flow, and we also expect the radiative efficiency to be close to 
the standard value for a NS, namely r] ~ 0.1 — 0.5. The spectrum of the radiation is also likely to 
be very different in the four models. This may provide a way to distinguish which if any of these 
possibilities is found in nature. 

In addition to the above possibilities, yet other flow conflgurations may be possible when 
we allow for the multi-dimensional nature of the flow. These could be explored with numerical 
hydrodynamics simulations. 

This work was supported in part by grants PHY 9507695 and AST 9820686 from the National 
Science Foundation. 
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Fig. 1. — Profiles of density p (g cm~^), proton temperature Tp (°K), electron temperature 
(°K), angular velocity 17 (in units of the Keplerian angular velocity at the NS radius Rns)} <^iid 
radial velocity v (in units of c) for accretion flows with a = 0.1 and (m, s) = (0.01,0.3) - solid line, 
(0.0001, 0.3) - short-dashed line, (0.01,0.1) - medium-dashed line, (0.01,0.7) - long-dashed line. 

Fig. 2. — The solid curve corresponds to the dimensionless specific angular momentum flux j, shown 
as a function of the dimensionless spin of the NS s, for m = 0.01 and a = 0.1. The long-dashcd 
line is j{s) for the thin disk case, taken from Popham &: Narayan (1991). The short-dashed curve 
corresponds to the analytical self-similar solution, equation (16). The lower panel gives a close-up 
of the small-s region of the upper plot. 

Fig. 3. — Variation of the angular momentum flux parameter —j as a function of m for s = 0.3, 
a = 0.1. The solid curve corresponds to the results from numerical computations, and the dashed 
curve corresponds to the analytical self-similar solution. 

Fig. 4. — Luminosity (erg/s) of the settling flow per unit logarithmic interval oi D = {R — 
Rns)/Rns for a = 0.1, s = 0.3. The curves correspond to different choices of the mass accretion 
rate: from below, m = 0.0001, 0.001, 0.003, 0.01, 0.02. A substantial part of the luminosity is 
independent of m; this part is from the settling flow. There is some weak dependence of the lumi- 
nosity on rh at small D, due to boundary layer emission, and at large D, due to non-self-similar 
behavior of the numerical solution. 

Fig. 5. — Same as in Fig. 1 for 7 = 4/3 and s = 0.3 (solid curve) and s = 0.01 (dotted curve). 
Corresponding self-similar slopes for an ADAF and the settling flow are also show. The long-dashed 
curves represent the same solution as the solid curve, but with ten times higher temperature at 
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